Constraints on decay plus oscillation solutions of the solar neutrino problem 
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We examine the constraints on non-radiative decay of neutrinos from the observations of solar 
neutrino experiments. The standard oscillation hypothesis among three neutrinos solves the solar 
and atmospheric neutrino problems. Decay of a massive neutrino mixed with the electron neutrino 
results in the depletion of the solar neutrino flux. We introduce neutrino decay in the oscillation 
hypothesis and demand that decay does not spoil the successful explanation of solar and atmospheric 
observations. We obtain a lower bound on the ratio of the lifetime over the mass of V2, (T2/m2) > 
22.7 (s/MeV) for the MSW solution of the solar neutrino problem and (T2/m2) > 27.8 (s/MeV) for 
the VO solution (at 99% C.L.). 
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, ^ , Solar neutrino experiments with Chlorine , Gallium ^ , and water cerenkov detectors [|| , , show unequivocally 
that there is a deficit of at the Earth compared to the predictions of the standard solar model |^ . It is commonly 
accepted that vacuum oscillations (VO) or matter induced MSW conversions, with large mixing angle (LMA) or small 
mixing angle (SMA), can account for the deficit of solar neutrinos. The available experimental results allow not only 
a test of the neutrino oscillation hypothesis but also offer possibilities of constraining new physics, e.g. neutrino 
i-H ' magnetic moment neutrino decay, flavour changing neutral currents, etc. Here we will be concerned with neutrino 
^ , decay in the context of neutrino oscillations. The possibility of solving the solar neutrino problem only through 
neutrino decay in vacuum was raised in Q and ruled out mainly owing to the fact that the lower energy pp neutrinos 
observed in Gallium experiments are less suppressed compared to higher energy ^Be and '^B neutrinos observed in 
Chlorine experiments. The suppression in neutrino flux caused by the solar matter induced decay to majoron has 
J correct energy dependence |^] but the required fast rates cannot easily be obtained in the standard scenario without 
OA ' conflicting with other constraints |^ . The oscillation plus vacuum decay scenario has been studied more recently in a 
two generation model in ^ , pO| . This analysis was prior to results from SNO Q . These results in combination with 
' earlier data from SuperKamioka ||^ can be used to separate the flux of the electron neutrino from that of other active 
Oh, flavours. This additional information provided by SNO can be quite useful in probing new physics. In the light of 
• this, we study in this paper the three generation model of neutrino oscillation plus decay including the recent SNO 
D , Q result in our analysis. 
i-C Radiative decays of neutrinos are severely constrained by laboratory experiments and a variety of astrophysical 

and cosmological observations (see references in the PDG [0). However, constraints on non-radiative decays are 
much less stringent. For a relatively heavy unstable neutrino, one can apply the argument that decay products would 
^ , contribute to the energy density of the universe and thus obtain a limit [|l2|. However, for masses of the order of eV 
■ or less, there is no limit on the lifetime. Solar neutrino data constrains non-radiative neutrino decay in this mass 
regime. 

We examine the scenario where the neutrinos from the Sun are depleted due to non-radiative decays like (a) Majoron 
emission decays 1/2 ^ i^i + J, where the i>i state is either sterile or the active antiparticle of vi, or (b) 1^2 ^ 3i^. We 
assume that the lowest mass neutrino i^i is stable, or at least that it has lifetime (ri/mi) ^ 20(s/MeV). (In the 
paper, we denote by rrii and the mass and rest frame lifetime of the Vi neutrino.) The neutrinos which arise as 
decay products in (a) or (b) can be either active or sterile. They are depleted in energy and carry on an average 
less than half of the solar neutrino energy. Since the flux of the electron neutrino sharply increases as we go down in 
energy, the relative contribution of the neutrinos produced in decay to the total signal in a given experiment will be 
quite small even if they are active. Because of this reason the net effect of the decay is depletion in the neutrino flux. 
This depletion is over and above the one caused by the neutrino oscillations or conversion which we assume to be the 
main cause of the solar deficit. The additional depletion caused by the decay is given by an exponential function 

where iJ^ refers to the energy of the decaying neutrino and R is the Sun-Earth distance. The solar neutrino survival 
probabilities are given in terms of the mixing si, mass difference Sm'^ = — m\ and the decay lifetime T2. We 
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determine the regions in the tan^ 0i — (5m^ plane allowed at 99% C.L. for different values of T2. The allowed region 
is found to shrink as T2 decreases and ultimately it disappears at some value which is taken to be the bound on the 
neutrino lifetime to invisible channel. The bound derived this way depends upon the specific solution of the solar 
neutrino problem and is stronger in case of the VO solution (r2/m2) > 27.8 (s/MeV) compared to the MSW-LMA 
and SMA solutions for which we get (T2/TO2) > 22.7 (s/MeV). In an earher pre-SNO analysis jl^, similar bounds 
were obtained for the LMA and vacuum solutions but there was no bound on T2 for the SMA solution. We show in 
the final section that inclusion of SNO result is critical for obtaining a bound on T2 . 

We assume the mixing among three neutrinos to be responsible for the deficit in the neutrino flux observed by 
the solar and atmospheric neutrino experiments. The mixing matrix U needed in order to accomplish this has the 
following form: 

= -C2SI C2C1 S2 ^2 (2) 

J V ^2Si -S2C1 C2 / yi^s J 

We have approximated U^a ~ in order to account for the negative results obtained in j/g disappearance experiment 
at CHOOZ In addition to the above mixing one also needs Sm'^ — — m\ « 10^'' — 10~^^ eV^ and 5a = 

'^3 ~ "^2 ~ 10~^ eV^ to account for the solar and atmospheric |]l^ neutrino scales. The required value of 5a is 
much larger than the value of the effective mass square 2^/2GFENe of the electron neutrino at the solar core. This 
suppresses mixing of the third neutrino in matter ||l^ which decouples from the rest in case of the MSW solution to 
the solar neutrino problem. The small value of Ue3 (taken here as zero) also prevents the mixing of the third neutrino 
in vacuum. As a result, the solar data cannot be used to constrain the lifetime of the heaviest mass eigenstate lying 
at the atmospheric scale. We therefore concentrate on limiting the 1^2 lifetime. 



II. VO PROBABILITIES 



Neutrinos are produced and detected as flavor eigenstates h'ajCt — e,fi,T but their time evolution operator is 
diagonal in the mass basis Vi^i = 1,2,3. The amplitude for flavor conversion during vacuum propagation is given by 



(3) 



where \Jai are elements of the mixing angle matrix and Ai (t) is the time evolution operator for the Vi mass eigenstate. 
For the decaying neutrino scenario the time evolution amplitude is given by 



t TJl 

Ai{t) = exp(-£'it) exp ( - — — ^ 



(4) 



where we allow for decay of the Vi mass eigenstate with lifetime Ti and we approximate Ei ~ p + mf/{2p). The 
probability for flavor conversion during propagation in vacuum at time t is then given by 



P^t) = Yj^c.^U*pU,pU*^^ cos 
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with El, = p. 

Now we assume that the lightest mass state vi does not decay. Using the mixing matrix 
conversion probability in vacuum from Sun to Earth is given by the expressions 
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where 
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= 1.5 X 10^^ cm — 500.3 s being the Earth-Sun distance. We notice that the decay hfetime does not occur in the 
expressions as Ue3 — 0. The charged and the neutral current interactions in the detector respectively are related 
to Pee and Pefj. + Per- Botli of these are seen to be independent of the atmospheric mixing angle. The sum of all three 
probabilities is also independent of the solar Sm^ , 

P=Pee+ Pep + Per = cj + sj exp ( - ^) (8) 



III. MSW PROBABILITIES 



If the flavor conversion of neutrinos in the Sun is by MSW mechanism then we have different expressions for 
probabilities for the decay plus conversion scenario. As discussed already, only two energy eigenstates Ui and 1/2 
corresponding to the lighter neutrinos participate in the MSW conversion. In the core of the Sun Ei > E2, there is a 
level crossing at the resonance point after which E2 > Ei i.e. ni2 > mi in vacuum. 

The probability of Ve vi just after level crossing is 

Pi = Pj-'^l + (1 - Pj)cL (9) 

Here, the first term stands for Ve going to V2 by mixing with probability and then jumping to vi with probability 
Pj at the level crossing. The second term means that Ve goes to vi by mixing and then does not jump to V2 with 
probability (1 — Pj) at the level crossing. The probability of Ve V2 just after level crossing is 

P2 = {1- Pi) = (1 - Pj)sl + {P.j)cl (10) 

In the formulas above, the Landau-Zener jump probability is given by 

^ expj-bsl/E^) - exp(-6/^^) 

b = ^ [ ) ~ lOM ^ ) MeV (12) 



\\A/A\resJ VeV- 



and A = 2^/2E^GFNe. The mixing angle in matter in the Sun is given in terms of the vacuum mixing angle by the 
expression 

...9« (-1 + 77(1-2^^)) 

= (l-2,(l-2.f)+ry^)V2 (13) 

^=^ = 6.6x10-^ A (14) 

After level crossing the vi state stays as it is but the 1^2 state can decay into antineutrinos or sterile neutrinos. 
Thus, at Earth the probability of detecting vi is Pi and of detecting V2 is P2 exp(— a/i?). We can now use the matrix 
(1) to find the Ve, J^fj,, v-r content of neutrinos at Earth, 

Pee = c^Pi + s\P2 exp(-a/£;,) 

= clsiPi + clclP2 exp(-a/S,) 
Per = slslPi + slclP2 exp(-a/^,) (15) 

As in case of the VO, both Pee and Pe^ + Per are independent of the atmospheric angle. But unlike VO, the sum of 
the three probabilities in (|l5| ) 

P=Pee+ Pet. + Per = Pi + (1 - Pi ) exp(-a/£;,) (16) 

depends on the solar scale through the Landau-Zener formula (pT|). 
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IV. EXPERIMENTAL RATES AND BOUND ON LIFETIME 



We use the probabilities derived above to obtain a bound on the neutrino hfctime. We consider only total rates 
for this purpose and include all experiments in our analysis. We do the analysis by two different methods. First we 
consider the total rates in Chlorine [Q , Gallium |^ , Super-K Q and SNO ^] experiments and determine the allowed 
regions in the tari^Oi — 5m? parameter space for different values of T2. The Cl^ Ga and the charged current rates of 
SNO by themselves cannot constrain T2 in the SMA region of the MSW solution since these experiments measure Pee 
which becomes independent of T2 in the small mixing angle limit. This is not true however in case of Pep, + Per which 
is probed by neutral current events in SK. Since the neutral current rates can be inferred by combining SK and SNO 
data , such combination is expected to constrain the lifetime T2. We show this explicitly by using the SK and SNO 
data alone. 

The rates of neutrino capture in the Chlorine and Gallium experiments can be written as 



J dE^^,(JaPe, 



i?„ = '-^^^ (17) 



where the subscript a = Ga, CI denotes the experiment and i = pp, ^Be, *B denotes the type of neutrino flux from 
the Sun. The spectra of the pp and neutrinos can be fitted with the analytical functions, 

= (5.95 X 10i°)[193.9(0.931 - £;^)((0.931 - E^f - Q.2&lY/^El] 
$B = (5.05 X 10^)[8.52 X 10^^(15.1 - E^f-'^^El] 

$Bc - (4.77 X IQ'')[5{E, - 0.862)] (18) 

where the neutrino fluxes arc in units of cm~^s~^ and Ey is in MeV. The first brackets in ( p^ ) give the total flux of 
neutrinos from the pp,Be, and B reactions and are taken from BP2000 [^, and the square brackets give the spectral 
shape I p^ . 

The Ga experiments can detect all three types of neutrino fluxes and the neutrino absorbtion cross section of Ga 
is given in |Q. The Chlorine experiment threshold is higher (0.8 MeV) and it detects only the Be and B neutrinos; 
we take the absorption cross section with the tables from ref. |l^ . 

The electron scattering reaction in Super-K and the charge-current deuterium dissociation reaction at SNO can be 
written as 

/ dE,a,^^nPee + J dE,G,^^n[P - Pee) 



and 



,CC _ / dE^OQQ^^Pe, 



The Vee" and v^.tC'' elastic scattering cross section after folding with the detector response function are tabulated 
in ref. ||l9|. The deuterium dissociation cross section are taken from pO |. 



Using the flux spectrum in equation ( |l8| ) and the cross sections |17|-20| we can calculate the theoretical rates for 
MSW or VO conversion probabilities as a function of the three unknown parameters: the 1/2 lifetime T2, 5vn? and 
the vacuum mixing angle 9i. The experimental rates R'^^^ with one-sigma combined (statistical and systematic) 
experimental errors Aq, are as follows 



R^^P^ = 0.335 ± 0.029 
R^^P^ = 0.584 ±0.039 
i?^^P* = 0.459 ± 0.017 
^SNO-cc = 0.347 ± 0.029 (21) 



From the theoretical Ra, and the experimental we compute the total for a-U experiments, defined as 
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( Z? Dcxpt \ 2 

X'^ E ^ ^ (22) 

a=Cl,Ga,SK,SNO " 



Setting the parameter a = we reproduce the standard contours of LMA , SMA and VO solutions shown as dotted 
contours in Fig. 1 and Fig. 2. For a = the global minimum of occurs in the vacuum region with x^j„ ^ 0.3. The 
Xmin + 11 contours which corresponds to 99% C.L. bounds for a are plotted for various values of non-zero a as shown 
in the solid contours in Fig. 1 and Fig. 2. We find that the LMA and SMA allowed parameter space disappear for 
a = 18 and a = 22 (in MeV units)respectively. The VO allowed region disappears for a = 18. These numbers 
translate into the bounds T2 > 22.7 s (r7i2/MeV) for the MSW solution (where we have taken the larger of the two 
bounds coming from SMA and LMA solutions) and T2 > 27.8 s (m2/MeV) for the VO solution. 

As discussed before the bound on T2 in the SMA region mainly comes from the SK data which probe Pe^ + Per ■ 
It is thus interesting to consider only the SK and SNO results which together determine Pee and Pe^ + Per in the 
same energy range (8 — 15 MeV). These two experiments are crucial in providing the bound on T2. To see the impact 
of these results , we plot the 1.96fT contours of the SNO (i?SNO = 0.347 ± 1.96 * 0.029 shown as dashed curves) and 
SK (i?sK = 0.459 ± 1.96 * 0.017 shown as continuous curves) rates. In Fig. 3 we show the allowed regions for a = 0. 
In Fig. 4 we plot the same SNO and SK rates for a = 25 and we can see that there is no overlap in the SMA region 
between the two experiments for this value of a. The bound on a obtained from SK and SNO is marginally weaker 
than the one obtained in the combined analysis of all (SK,SNO, Ga and CI) experiments. 

In summary, we have used available results of the solar and atmospheric neutrino experiments to constrain neutrino 
lifetime within the three generation picture of neutrino oscillations. While the exact bound depends upon the specific 
solution, we typically find ^ > 28 {-ff^) The corresponding bound T2 > 2.8 • 10~^sec for eV mass neutrinos is not 
very strong but useful since it is the only one following from laboratory experiments in this mass range. 
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FIG. 1. Allowed parameter space for VO plus decay solution, continuous curve is for decay parameter a = and dashed 
curve is for a = 10. 
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FIG. 2. Allowed parameter space for the MSW plus decay solution, continuous curve is for decay parameter a = and 
dashed curve is for a = 10. 
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FIG. 3. The 1.96(j allowed parameter space for Rsk as region enclosed by dashed curves and Rsno as region enclosed by 
continuous curves; for decay parameter q = 0. 
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FIG. 4. The 1.96(7 allowed parameter space for Rsk a.s region enclosed by dashed curves and Rsno as region enclosed by 
continuous curves; for decay parameter a = 25. There is no overalap in the SMA region for a = 25. 



